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Abstract 



In this paper we study a novel means of coupling neutrinos to a Lorentz violating 
background k-essence field. We first look into the effect that k-essence has on the 
neutrino dispersion relation and derive a general formula for the neutrino velocity in the 
presence on a k-essence background. The influence of k-essence coupling on neutrino 
oscillations is then considered. It is found that a non-diagonal k-essence coupling 
leads to an oscillation length that goes like A ~ E^ 1 where E is the energy. This 
should be contrasted with the A ~ E dependence seen in the standard mass-induced 
mechanism of neutrino oscillations. While such a scenario is not favored experimentally, 
it places constraints on the interactions of the neutrino with a cosmological k-essence 
scalar background by requiring it to be flavor diagonal. All non-trivial physical effects 
discussed here require the speed of sound to be different from the speed of light and 
hence are primarily a consequence of Lorentz violation. 
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1 Introduction 



The observation of the accelerated expansion of the universe [1, 2] has been one of the most 
important recent discoveries in cosmology. Many possible explanations have been put forward 
which may be classified under two general classes: models with a cosmological constant, or 
dynamical models of dark energy. Their common feature is to provide fluids with negative 
pressure to drive the acceleration. Among the dynamical dark energy models, only k-essence 
[3] , has the advantage of explaining not only the current phase of accelerated expansion but 
also the coincidence problem, i.e., why the cross-over from the matter dominated era to the 
current era happened so recently in the past. This explanation, however, is not without 
its own problems as was first pointed out in [4]. It was shown there that in order to solve 
the coincidence problem, the universe had to go through an era where the sound speed or 
the speed of propagation of the k-essence fluctuations must become super-luminal. This 
problem was addressed in [5], where it was shown that at the classical level, super-luminal 
propagation does not necessarily imply causal paradoxes. In particular, propagation in a 
k-essence background does not have any additional causal difficulties over general relativity 
where the only problems are associated with space-times which admit closed time-like curves. 
In the course of this analysis, the authors found a very interesting way of describing the 
propagation of k-essence fluctuations in terms of an emergent metric which depends not 
only on the space-time metric but also the background cosmological k-essence scalar. One 
may thus think of this nontrivial, Lorentz violating cosmological k-essence background as 
the aether in which matter perturbations propagate. This emergent metric description is 
used in this paper to couple the k-essence background to neutrinos. 

If there is a scalar field pervading the universe then the effective field theory viewpoint im- 
plies that it must undergo interactions with the matter that is present. The question of the 
observabilty of dark energy directly through its couplings to ordinary matter is an important 
one [6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17] and this paper attempts to address aspects of 
it in the context of k-essence. On short distance scales the universe is inhomogenous with, 
in particular, plenty of black holes. Thus the interactions with dark energy with black holes 
could be one way to study the above mentioned question. In this paper, however, we will be 
concerned only with the effect of the k-essence background aether on the propagation of neu- 
trinos. The coupling of fermions to this background is in itself an interesting question from 
the theoretical point of view. Most studies of the interaction of dark energy with fermions, 
couple them through a Yukawa like interaction [16], which is quite reasonable. However, in 
this paper, we do this differently using the vierbeins constructed out of the emergent metric. 
Throughout this article the k-essence field is treated strictly as a background, however we see 
no reason not to treat it as a dynamical field. As we argue in section §2, within the effective 
field theory methodology, the terms with higher derivatives of the fermion fields do not give 
rise to ghosts. The main focus of this paper is on looking for observable consequences of dark 
energy, so in this case the fermion in question is the neutrino which we show undergoes flavor 
oscillations when traveling through the k-essence aether. However, the way we introduce the 
fermion-k-essence coupling could be used to obtain new types of interactions between dark 
energy and other forms of matter including dark matter. The emergent metric from [5] which 
is used throughout this paper to couple the dark energy to neutrinos is covariantly constant. 
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In a future publication we will discuss how to consider even more generalized couplings by 
introducing torsion in the emergent space-time. 

It has been noticed previously [5], in the context of the propagation of k-essence fluctuations 
in a classical background that defines the emergent space-time, that Lorentz invariance is lost 
when the speed of propagation of the fluctuations (the speed of sound c s ), is different from the 
speed of light. The same is true for the propagation of neutrinos in this background. In fact, 
all of the physically interesting results that we obtain in this paper are present only when 
c s 7^ c, i.e., when there is Lorentz violation. Nontrivial neutrino flavor oscillations require, 
in addition, non-diagonal flavor couplings of neutrinos to the k-essence background. In the 
past various models of neutrino oscillations have been considered which require an explicit 
violation of Lorentz invariance [18] or of the equivalence principle, [19, 20]. The energy 
dependence of the oscillation length is the same in these models as the one considered in this 
paper. In this sense our model may be considered a theoretically and phenomenologically 
motivated manifestation of the same phenomenon. We should emphasize that we do not 
have any violation of the equivalence principle, but the emergent metric which contains 
contributions from the k-essence background can be different for different flavors of neutrinos. 
This is made possible in our model by a flavor non-diagonal coupling in the part of the 
emergent metric involving the k-essence background only. 

K-essence is a theory with noncanonical kinetic terms and coupling it to neutrinos through 
the vierbein of the emergent metric alters the speed of propagation of the neutrinos. The 
consequent dispersion relations are analyzed in section §2. The data from supernova 1987 is 
then used as an input to constrain some of the parameters of our model. In section §3 we 
consider the possibility of neutrino oscillations induced by their coupling to the k-essence 
background. The more interesting and novel case is when the different neutrino flavors couple 
with different strengths to the k-essence field. Oscillations are induced essentially due to the 
fact that the speeds of propagation of the different neutrino species in the background aether 
are consequently different. In this case neutrinos would oscillate even if they were massless. 
In section §3.2 we discuss this case in some detail and obtain a general formula for the 
oscillation probability with massive neutrinos. Our results are quantitatively different than 
the case of flavor oscillations with only massive neutrinos. In particular, the oscillation 
length varies with the inverse power of the neutrino energy. Such a behavior is ruled out 
by the data from the Kamiokande experiment [21]. Thus, we are able to place bounds on 
the allowed strengths of the k-essence coupling to neutrinos. In particular, the data strongly 
favor diagonal flavor couplings of k-essence to neutrinos. As a preliminary to this analysis, 
in section §3.1 we discuss the case of neutrino oscillations with massive neutrinos but with 
equal coupling strengths of all flavors to the k-essence background. Here we find a rather 
simple modification of the well known formula for the flavor oscillations of massive neutrinos 
with the only difference arising due to the fact that the neutrinos travel along geodesies in 
the emergent space-time. In section §4 we present our conclusions. Certain technical details 
of the coupling of Dirac fermions to the emergent metric are relegated to an appendix. 



3 



2 Neutrino Coupling To a K-essence Background 



Before we can talk about neutrino interactions with k-essence we should first review the 
latter [3, 5]. In general k-essence is a theory of a scalar field with non-canonical kinetic 
terms. The lagrangian of a single k-essence scalar is usually denoted by a single function 
L(X, 0) where X = \ V^V^. For a given solution to the k-essence equations of motion, 
the behavior of perturbations \ = 8<f> in the k-essence field around the background can be 
written can be described by a canonical scalar field action, but with the space-time metric 
g^u replaced by an emergent space-time metric G^ given by [5], 

G, v = JT 2 (g^ + ^I^V^V^ (1) 
where Q 2 = jf^ 1 and c s is the sound speed of k-essence fluctuations,which is defined as 

Here p& and p^ denote the energy density and pressure of the k-essence background, and the 
subscript signals that (2) should be evaluated while holding constant. In the case where 
V^0 is a time-like vector (i.e. X > 0) then it can be shown that pk = L and pk = 2XLx — L, 
and therefore it follows that 

s L x + 2XL X x ^ 

Because g^ is replaced by in the \ field action, the characteristics of x follow the 
geodesies of G^ and not g^ u . This interesting fact implies that x has a different causal 
structure then all other fields. In particular, if in some frame g^ u = r]^ and the background 
k-essence field is uniform, then the emergent space-time metric is given by 

G^dx^dx" oc c 2 s dt 2 - dx 2 (4) 

The metric G^ u defines a different causal structure than g^ u . The "light" cones of G^ are 
defined by characteristics with velocities c s instead of the speed of light. It was shown in 
[5] that even if c s exceeded the speed of light, causality in a k-essence theory would still be 
preserved despite the superluminal speed of k-essence fluctuations. This fact can be roughly 
understood by thinking of the k-essence metric (4) as the space-time interval in the special 
relativity but with a different value for the speed of light. Therefore causality is preserved 
in k-essence for much the same reason that it is preserved in special relativity. However, in 
order for this rational to hold the k-essence action most satisfy certain constraints. These 
constraints arise from the need to make G^ have the proper signature. This translates to 
the requirement that the k-essence action satisfies 

l + ^^>0. (5) 
Lx 



1 Derivatives of L with respect to X are denoted by a subscript, so that Lx = §jf and Lxx — -qx? 
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Note that this condition is equivalent to the stability constraint: c 2 s > 0. 

In this paper we will take inspiration from the k-essence perturbation action, and consider 
the possibility of other fields coupling to G^ v . In particular, we will take the action of a 
neutrino coupled to a gravitational metric g^ v and replace this with the k-essence metric 
Gpy. The action of our hypothetical k-essence-coupled (Dirac) neutrino v is given by: 

S = J d^xEv [i^V„ - M] v (6) 

where E = detE^ and 7 M = E£ / -f a , and 7° are the standard gamma matrices. The vierbein 
field E% of the emergent space-time geometry and its inverse are given by 




Here we have defined u = j- — 1 and u = —j^^ 2 - Further note that = g^ u V u (j). Also, 
we have included a coupling constant g, which accounts for the interaction strength of the 
k-essence background to the neutrinos. The emergent metric G^ v is still covariantly constant 
as in [5]. 

At this point we would like to emphasize two features of the action. First, we note that when 
c s ^ 1 the model is not invariant under Lorentz transformations [5]. As we will see in detail 
in section §3, all the physical effects that we discuss in this paper are consequences of this 
Lorentz violation in the sense that they vanish at c s = 1. Other features like nondiagonal 
flavor couplings are also important to get nontrivial flavor oscillations, however, Lorentz 
violations must always be present. Secondly, in this model we treat the k-essence background 
as a classical field, which does not experience any appreciable back reaction effects from the 
neutrino field. This allows us to treat the neutrino/k-essence coupling as a contribution to 
the kinetic term in the neutrino action. If the k-essence field was dynamical, this would 
lead to higher order derivatives of the neutrino field in the k-essence equation of motion, 
which could potentially create ghosts in the quantum theory. This is not a problem for 
this paper since we treat the k-essence scalar strictly as a classical background, however, 
we would like to emphasize that within the effective field theory methodology treating the 
k-essence scalar as dynamical would not give rise to such problems in any case. The terms 
with higher derivatives of the fermion field would be considered as higher order in the low 
energy effective action expansion. 

The derivative operator P M in (6) represents the spinor covariant derivative with respect 
to the emergent k-essence background described by E%. The proper definition of X> M and 
its specific form in the case of a general k-essence field in a flat space-time background are 
given in appendix A. From here on out we will ignore the spinor connection term in "D M , 
which we justify on the basis that higher derivatives of the k-essence field are negligibly 
small at the present time in most models of k-essence. However, if at one time in the history 

2 Note, that in our definition of u and u we have assumed that the 4-gradient of the k-essence background 
field is time-like. In the case of space-like V^, the sound speed is given by the inverse of the expression for 
the sound speed in the time-like case. Thus, for space-like one should replace c s — > — in the definition 
of u and u. 
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of the universe the connection was not small, its effect can be accounted for by replacing 
M 2 — > M 2 + R/4, where R is the scalar curvature of vierbein E^. The emergent space-time 
metric is by definition given by G^ v = E a E h v r\ ab , which is 



= E«E b uVab = n~ 2 [g, u + V M 0V^ 



where w = 2gu + g 2 u 2 . Note that if g = 1 then w — c 2 — 1, which is the standard result for 
the emergent metric in k-essence theories. The inverse of G^ v is given by 

G» v = Q 2 (g^ + ^ V ^ V! >) ( 9 ) 
where w = 2gu + g 2 u 2 = — y^=. The determinant of G^ v is given by 

-(det^)^ (10) 



E 2 = - det G^ 



In order for this modified k-essence induced metric to have the proper signature: 1 + w > 0. 
Note that this condition is equivalent to (5) in the case when g — 1. The k-essence coupled 
Dirac equation reads 

{i^V^-M)v(t,x) = (11) 
If we square this equation we can obtain a Klein-gordon equation for z/(t, x) 

{G^VpVv + R/4 + M 2 )is(t, x) = (12) 
Here, R is the curvature scalar whose definition in terms of the connection form is 

r = [d,Q ab u - d v w\ + w c ^\ - n%,n cb li ] (is) 

It is evident from (12) that the curvature scalar acts as a mass term. However, as we 
discussed earlier, in most k-essence models higher derivatives of the field <fi will be negligible, 
and therefore we can ignore R from here on out. By taking the plane wave approximation 
for the neutrino field the phase of v is proportional to e - ' l J rp ^ dxh ' . If we assume that the 
interaction of the neutrino field is weak and the background geometry is flat, then the 
dominant space-time dependence of the neutrino field comes from the phase factor. Thus 
the Klein-Gordon equation in momentum-space leads to the dispersion relation 

G^ PflPu — M 2 = (14) 

Define an effective momentum p M as 

p, = n-^Efa, = p, + H(^V»V M (15) 

The covariant and contravariant effective momentum are defined with respect to the space- 
time metric and not the emergent k-essence metric G^. Thus the index on p^ is raised and 
lowered using the space-time metric g^ u . Because of this property, it follows that 

n%p» = Q 2 g^p,p v = g^(e«Ee Pp )(e b u E£px) = Tf*E*E% Pl n, = CTp^ (16) 



6 



For the purposes of this paper we will assume that the background space time is flat so 
that g^ v = r]^ u and = 5%. This is in fact, a good approximation cosmologically for the 
applications we have in mind in this paper. Therefore, the dispersion relation in terms of the 
effective momentum is 

PoP°+PiP l = m 2 (17) 

where m 2 = Q~ 2 M 2 is an effective mass that we have defined here for convenience. We wish 
to use the on mass shell condition above to find the particle velocity v of the neutrinos which 
is represented by the group velocity v = J^. In general we will find that the velocity of 
a neutrino depends on the energy which stands in contrast to the case of a free massless 
neutrino, where the velocity is energy independent. 

Throughout this paper we will define p^ = (E, — p), p = |p|, p = ph, and = Vo0. In the 
next few sections we will find the neutrino velocity for two special uniform k-essence 

field, and a static k-essence field, after which we will derive the velocity assuming the most 
general k-essence field configuration. 



2.1 Simple Case: is Uniform 

Before we try and find the expression for the neutrino velocity with the most general k-essence 
field, let's find the velocity when is uniform. This is probably the most relevant case since 
most k-essence theories, in particular those that attempt to address the cosmological constant 
problem, assume that the spatial derivatives of the k-essence field are negligible compared 
to its time derivative [4, 3, 22, 23]. If is uniform then Vj0 = 0, and thus 



pop + pip 1 = m 2 =>- E - 



1 + 



gu(j> 2 

1x~ 



By definition X = \<fi 2 - Therefore 




(19) 



and the group velocity, which represents the neutrino particle velocity is 



where c v = = nz§t +g - The speed c v plays the same role for massless neutrinos that 
the sound speed c s does for massless k-essence perturbations. Both c u , and c s represent 
limiting speeds that neutrinos and k-essence perturbations, respectively, are required to not 
exceed as measured from the frame in which the k-essence background is uniform. The 
fact that in general c v ^ c s is due to our inclusion of an arbitrary coupling parameter g. 
Phenomenologically this coupling is small, so we easily see that if c s > 1, then so is c v . 
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If m 2 = then the neutrino velocity (20) will be equal to c v . As a massless particle coupled 
to the emergent background geometry, the neutrino will travel on the null geodesies of G^ u 
not g^ w . In the uniform case, the emergent metric (8) with arbitrary g is 

G^dx^dx" oc c 2 u dt 2 - dx 2 (21) 

As one can see here, null lines in travel at a speed c v . In this sense the emergent geometry 
of a uniform k-essence field acts just as a minkowski space-time except that limiting speed 
is now c v instead of the speed of light. 



2.2 Slightly Less Simple Case: is Static 

In direct contrast to the last case let's consider what happens to the neutrino velocity when 
<fi is time independent, but has non-zero spatial gradients. In this case V M is a space-like 
vector (i.e. X < 0) and the energy density and pressure are instead given by p = —L and 
p = L — 2XLx- It follows that c 2 = Lx+ ^ xx , which means that the definition of the sound 
speed for a space-like k-essence field is the inverse of sound speed for a time-like k-essence 
field. Thus we make the replacement c s — > — . It follows therefore that u = c s — 1 now. 

C-s 

We should note, however, that the metric has not changed, only the definition of c s . If we 
assume a static but spatially varying k-essence field, then the on mass shell condition states 
that 

PoP° + PiP i = m 2 E 2 = p 2 [l + w cos 2 9] + m 2 (22) 

where 9 is the angle between p and V0. Note that in this case X = — ||V0| 2 . The neutrino 
phase velocity is then given by 




and the group velocity is 



(23) 



1 + w cos 2 9 . . 

v 9 = i — ( 24 ) 

/ l + wcos 2 ^ + ^ 



V 



It is informative to evaluate (24) at the two extremes of cos 2 9; that is when p and V0 are 
parallel and when they are perpendicular. If p and V0 are parallel then the angle 9 between 
them vanishes and we find that the neutrino velocity is 

v(B = 0)= , 1+9U „=^=^ (25) 



where now c v — 1 + gu — c s g — g + 1. Again, if g — 1 then c u = c s . We can see that 
the velocity of the neutrino (25) is almost the same as the formula given in (20) except 
that the mass term in the denominator now has a factor of \. Again if the effective mass 
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of the neutrino is zero then the neutrino velocity is equal to the sound speed of k-essence 
fluctuations. This is due to the fact that the neutrino propagates on null geodesies in the 
emergent k-essence background. Now, on the other hand, if p and V0 are perpendicular 
then 



V {Q = tt/2) = 




and the formula for the velocity of the neutrino is the same as it would be if in the absence 
of a k-essence field. This is because in the static field case the coupling between the neutrino 
and k-essence is proportional to p • V0. This means that the neutrino will act as a free 
particle propagating on a flat lorentzian space-time whenever it is traveling perpendicular to 
the direction of the field gradient. The directional dependence of the neutrino velocity in a 
spatially varying k-essence field stands in stark contrast to the neutrino velocity in uniform 
k-essence. If neutrinos traveling from a distant galaxy were to travel through a region of 
spatially varying but static k-essence field (a k-essence halo) on their way to a detector on 
earth, we should expect to see evidence of an-isotropy. However, even if it were possible to 
detect a sufficiently large neutrino flux, it is expected that any spatial variation of the k- 
essence field will be very small compared with its variation in time. Therefore any anisotropy 
in the neutrino velocity would most likely be unobservable. 



2.3 Neutrino Velocity In a General K-essence Background 

Without making any assumptions about the nature of the k-essence field, the on mass shell 
condition (17) becomes 



w<fi 2 E 2 w • — ,.E w _ ,. 2 ^ m 2 



The solution for the — is 

v 



E_ ~#0(n ■ V0) ±y/l + #(0 2 - (n ■ V0) 2 ) + (1 + g)g 
p 1 + ^ 



(28) 



The choice of either plus or minus sign in the solution reflects the two particle/anti-particle 
states of the neutrino, with the plus sign corresponding to the neutrino and the minus sign 
corresponding to the anti-neutrino. According to the Feynman-Stueckelberg interpretation 
of anti-particles, the anti-neutrino can be thought of as a positive energy neutrino traveling 
backwards in time. Thus, the solution (28) with the negative sign, representing the anti- 
neutrino energy, should have a overall negative sign removed. Furthermore, since time is 
reversed, this means that in order to have the anti-neutrino traveling in the direction of 
n, we must replace n — > — n when we choose the negative sign in (28). In the end, the 
energy-momentum relation for the neutrino and anti-neutrino will be the same and given 
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by (28) with the positive sign. If we expand (28) to first order in u and m 2 , then for the 
(anti-)neutrino becomes 

E _ gu f • A _ , \ 2 m 2 
2X 

From (28), we find that the group velocity is 



P + +v (29) 




1 + ^ \ «/1 -I- -2LU2 

It is important to note that the neutrino velocity (30) does not change under a redefinition 
of the k-essence field variable <fi unless m 2 ^ 0. It is easy to show that under a redefinition 
from <fi to another field tp defined by = g((p), then (30) would remain unchanged were it 
not for the the m 2 term in the denominator. Recall then m 2 is not the physical neutrino 
mass but rather a rescaled mass, which is rescaled by the conformal factor Q in the emergent 
metric. After a field rescaling the effective neutrino mass becomes 

9 M 2 ~L X o M 2 Lj> m 2 

fi2 Cs [«^)]2 Cs [^(^)]2 ^ ^ 

Before any objections are raised we should point out that from the beginning we have chosen 
a specific background k-essence field that the neutrino couples to. In essence what we have 
done is to fix the "gauge" of the k-essence field. Therefore, it is no surprise that by changing 
the field variable we are changing the physics of the neutrino field. Since a field redefinition 
changes that conformal factor Q, a field redefinition is itself a conformal transformation. That 
the effective mass is the only quantity that changes under a field redefinition is a reflection of 
the fact that by adding a neutrino mass we are in essence breaking the conformal invariance 
of the neutrino action. 



2.4 Comparisons with Observation 

In 1987 a supernova was observed [24, 25] in the Large Magellanic Cloud that provided a 
limited but unique opportunity for the study of neutrino physics. During this event three 
neutrinos where detected at neutrino detectors here on earth and unambiguous identified as 
having been the emitted by the supernova. By comparing the time interval between when 
the supernova was first seen and when the neutrinos were detected a limit on the neutrino 
speed's deviation from the speed of light could be bounded. In [26] the authors calculated, 
using the available data from the supernova event, that that deviation of the neutrino speed 
from the speed of light can not be more than 1 part in 10 8 . In order words if v v is the 
neutrino speed then 
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< 10~ 8 . (32) 



If the neutrino is massless then v u = c u , and the left hand side of (32) is equal to \gu\ in 
the physically relevant static case 3 . Thus, in the massless neutrino limit (32) represents an 



3 note that c = 1 in all previous sections 
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observationally required upper bound on \gu\. In the most general case this bound becomes 
(setting c = 1 again) 

< 10~ 8 (33) 

The most generous upper bound that can be placed on gu is \gu\ < 10~ 8 . This can also be 
translated into a restriction on the k-essence sound speed c s . If the k-essence field has a time 
like gradient (i.e. X > 0) then the range of values for c s is 



i + io-VbP i - io-«/bl 

As of this moment there is an insufficient amount of data to constrain c s from cosmological 
observables. Recent studies have shown that fits of general dark energy models to the current 
CMBR data are largely insensitive to the value of c s [27, 28]. Without any other observational 
constraints on (or better yet a value for) c s , it is impossible to reliably estimate the value of 
the k-essence coupling g. However, with more precise data in the future it may be possible 
to get a better handle on the value of c s , and once that has been established, (34) can be 
used to put useful restrictions on g. 
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< 10" 



3 Neutrino Oscillations 

Experiments [21] have confirmed the the phenomenon of flavor oscillations whereby neu- 
trinos oscillate between various possible flavor eigenstates as they travel away from their 
source. There are different ways of explaining this oscillation, but all mechanisms for induc- 
ing neutrino oscillations involve some term in the neutrino lagrangian that is non-diagonal 
in the flavor eigenstates. Although the most popular way for inducing neutrino oscillations 
is by introducing a mass term [29], several other mechanism have been proposed over the 
years, such as: violation of the equivalence principle (VEP) [19, 20], torsion induced neutrino 
oscillations [30], violation of Lorentz invariance (VLI) [18], violation of CPT symmetry [31]. 

If neutrinos do indeed couple to a k-essence background in the manner we described in 
§2, then it is possible that k-essence can play a role in neutrino oscillations. There are 
two ways in which k-essence could effect neutrino oscillations. If the k-essence coupling is 
the same for each neutrino flavor then the energy difference between energy eigenstates will 
not be affected. However, neutrinos coupled to the k-essence background will travel along 
geodesies in the emergent space-time. This will have an affect on the phase of the neutrino 
wavefunction, which will be observable in the neutrino oscillation probability [32, 33, 34, 35]. 

Another way k-essence can influence neutrino oscillations is if the k-essence coupling g is 
non-diagonal in the flavor eigenstate basis. Imagine now a model of two neutrino flavors that 
couple non-diagonally to k-essence in the flavor eigenbasis. The lagrangian of this 2-neutrino 
system can be written as (a, (3 = 1,2 are flavor indices), 

£ = i ^2 E a U Ka j£V fjPKa ~ ^ VKa ( Ma P E P + E a M afi) VRf3 (35) 

a=l,2 a,P=l,2 
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where E a = E aa and E a/3 = [det E^(g)] a p. In this lagrangian the k-essence coupling g has 
been replaced by a matrix valued object g, that is not necessarily diagonal in the flavor 
and mass eigenstates. We have defined v K as the "k-essence eigenstates" , which are the 
eigenstates of the k-essence coupling matrix g. In general the k-essence eigenstates will not 
be the same as the neutrino flavor eigenstates. Because of this non-diagonal coupling of 
k-essence to the flavor eigenstates, the formula for the oscillation probability in the case 
of k-essence induced neutrino oscillations (KINO) will differ significantly from the typical 
mass- induced result. 

In this section we will assume that the k-essence field is weakly varying, with small second 
derivatives so that we may effectively treat the k-essence interaction, which only involves the 
first derivatives of the scalar field, as a coupling constant. This greatly simplifies matters, 
because it allows us to diagonalize the neutrino equations of motion in the momentum space 
representation. This assumption is consistent with the literature where most models of 
k-essence take the field and its sound speed to be relatively constant in time and space. 
As a result of this assumption we can absorb the determinants E a into a redefinition of 
the neutrino wavefunction vk and mass matrix M a p. Therefore, it is safe to ignore this 
contribution in our analysis. 

In the first part of this section we will consider the effect of k-essence on neutrino oscillations 
when the k-essence coupling is equal for each flavor eigenstate. 

3.1 Neutrino Oscillations with Flavor Diagonal K-essence Cou- 
plings 

At first it may seem that when both neutrino flavors couple to k-essence identically, the usual 
formula for mass -induced neutrino oscillations should not change. However, even in this case 
the fact that neutrinos travel on geodesies in the emergent space-time will create deviations 
from the flat space result. In order to best analyze neutrino oscillations in the presence of 
a K-essence background, we will use the simple method used in [32, 36] to study neutrino 
oscillations in curved space-times. In all cases, the important quantity of interest when 
calculating the neutrino oscillation probability is the phase of the neutrino wavefunction. 

\u f {t,x)) = e~* \u f ) = e-'S^ \ Uf ) (36) 

The difference between the phases $ of two different flavors of neutrino will lead to neutrino 
oscillations. In the expression (36), the hats over <E> and the 4- momentum indicate that 
these are operators which act in the flavor space of neutrinos. The phase $ when written 
out in terms of the energy and three momentum is 

$ = / Eldt - p • dx (37) 
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where I is the identity operator 4 . In a full and proper treatment of neutrino oscillations, the 
neutrinos must be modeled by spatially localized wave packets composed of neutrino energy 
eigenstates. However, if one uses this approach to study neutrino oscillations, they would 
find that in the relativistic limit it is acceptable to assume that both neutrino eigenstates 
propagate on the same null geodesic between the neutrino emitter and detector in space-time 
[37, 38, 39] 5 . Therefore the phase operator (37) becomes 



$ = / [ EI — ft 



~dt 



dt (38) 



where t e and td are the values of the coordinate time at which the neutrino signal is emitted 
and detected. The "0" subscript on ^ denotes that this quantity is to be evaluated along a 
null geodesic between the emitter and detector. To find we start from the definition of 
the canonical momentum of a massive neutrino 

dx v 

Pfl = mG^— (39) 

Here s is a proper time coordinate defined in the neutrino rest frame. With this, it is clear 
that ^ is 

at 



rfx _ g _ g% P - %(<}>E + p ■ V0)V0 
dt | G 0v p u E + £(j>E + p-V<l>)j> 



(40) 



Along a null geodesic the energy and momentum of a massless neutrino are denoted by E^ 
and Let us denote the unit vector in the direction of the neutrino momentum by n. 
Thus, n • [^] is given by 



n 



~dt 



r (fi.v^) + vi + s(^-(A-vOT E m 



J 



-i _|_ w <f> 2 p(o) 

1 ^ 2X 



To relate E and p, recall the discussion from section §2.3. It is easy to show from the work 
done there that the momentum as a function of E expressed to leading order in m 2 is 



^(n-V0) + A /l + ^(^-(n.V0) 2 ) 

p w E i w (h / — = (42) 

l- 2X (n V0) / 1+ -( 2_ (fl . V 0)2) 



2X ' 

Now remember that m 2 = Q~ 2 M 2 . In order for there to be neutrino oscillations the flavor 
and mass eigenstates must not be the same. Therefore, we have to replace the M 2 with a 



Note, we could have assumed that the momentum operator was proportional to the identity and the 
energy E was an off diagonal operator. Both approaches are equivalent when a first order expansion in the 
energy has been taken, as it will be in the final answer. 

J In the wavepacket treatment, an exponential damping is found in the final result for the oscillation prob- 
ability in the relativistic limit. This is due to the decoherence of the superposition of neutrino wavefunctions. 
So far there has been no solid evidence for decoherence effects [40] , and we are therefore safe in ignoring this 
possibility in other analysis 
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matrix M 2 that is non-diagonal in the flavor eigenstate basis. Since only phase differences 
between flavor eigenstates will be important, we can ignore terms not proportional to M 2 
in the phase operator (38). With (41) and (42) the expression for the phase operator (38), 
modulo terms proportional to the identity, becomes 



M 2 L( 



2E<® 



(43) 



where x e and Xd are the neutrino and detector positions, respectively. We have defined an 



effective distance L, which is a function of the coordinate distance I 
neutrino emitter and detector. The expression for L(£) is given by 



I %e 



L(£) 



■f t -/ 



n- 2 (h-dx) 



l + ^(0 2 -(n-V0) 2 ) 



between the 



(44) 



In general this distance will differ from the true emitter/detector separation due to the 
neutrino coupling to the k-essence medium. There are two sets of neutrino eigenstates in 
this system: the flavor eigenstates |z//) = {\u e ) , (z/^)} which couple diagonally to the weak 
current, and the mass eigenstates \v m ) = {|^ mi ) , |^m 2 )} that define the basis in which the 
mass matrix is diagonal. In general these two sets will not be equivalent. However, we can 
define an SU(2) matrix V such that 



\ U f) =^2 V fmWm) 



(45) 



The matrix V is referred to as the mass mixing matrix. The most general SU(2) matrix can 
be represented in terms of an angle and three phases: 











cos 9 sin 9 
— sin 9 cos 9 







(46) 



In general, V is an SU(2) matrix and will therefore have three additional phase degrees of 
freedom. However, in this case we can absorb these phases into a redefinition of the neutrino 
wavefunction. Thus we can safely ignore these phases in V and as a result the mixing between 
the flavor and mass eigenstates is determined by a single angle which we will call 9m- 



V 



cos 9 m 
— sin 9 m 



sin 9 M 

cos9 M 



(47) 



Written in the flavor eigenstate basis the phase operator (43) is 



2E 



Am 2 L{ 



- cos2#a/ 
sin 29m 



cos 26 m 



+ 



m 



2E 



1 
1 



(4£ 



where mi and m,2 are the masses of the two neutrino mass eigenstates, and Am' 



VTLry 



and m 2 



Note that we have dropped the "(0)" superscript on the energy E merely 



out of convenience. Since we are interested in finding probabilities then we can subtract 



14 



from the phase matrix any term proportional to the identity matrix without changing our 
final result. Therefore, we will ignore the very last term in (48) from here on out, and the 
phase factor in (36) becomes 



-i sin 



,Am 2 L(£). 
~AE 



— cos 29 m sin 29 u 
sin 29m cos2#m 



,Am 2 L{ 



cos 



AE 



1 
1 



(49) 



The factor e -1 * plays the role of a evolution operator. Denote by \vf) the initial state of a 
neutrino of flavor /. The time at which the neutrino is created can be defined, without loss 
of generality, as occurring at t = 0. The wavefunction of this neutrino at some later time 
t > after it has traveled a distance £ away from its source will be given by the action of 
e -»*W on the neutrino wavefunction \vf): 



M«)> 



/' 



(50) 



Applying the phase operator (49) to the neutrino ket \ ff), we find that the flavor eigenstate 
wavefunctions at a later time t > when the neutrinos have traveled a distance £ are 



M4> = - 



,Am 2 L(£) w nn 
i sin l — ) (cos 29 M \ u e ) 



-i sin 



AE 
,Am 2 L(£) 

~AE 



sm26' M + cos( — — — ) |i/ e ) 

4£< 



Am 2 L(£) 

(sin 20 A/ |i/ e ) + cos 26» M |^)) + cos( — — ) |^) 

4£/ 



(51) 
(52) 



Suppose a neutrino of flavor / is created from some source of interest (e.g. the sun, atmo- 
sphere, nuclear reactor, etc..) at some time t e . The probability of this neutrino appearing 
as an /' flavored neutrino in a detector here on earth at some later time after having 
traversed a distance £ is given by the expression 



P(u f -> v r ) = | {v f \ u f (£)) | 2 = 5 fJ i + (-1) 5 /./' sin 2 20 M sin 2 ( 



Am 2 L(£), 
AE 



(53) 



As we can see here the formula for the oscillation probability is almost exactly the same as 
in the standard mass-induced result. The only difference is the effective distance L is used 
instead of the coordinate distance £. 



3.2 Neutrino Oscillations with Flavor Non-diagonal K-essence Cou- 
plings 

At the beginning on this paper we have introduced a k-essence/neutrino coupling parameter 
that we have denoted by g. In the last subsection we assumed that this coupling was the 
same for each neutrino flavor. In this subsection, however, we will consider what happens 
when this coupling g becomes a matrix valued object g that operates in the neutrino flavor 
vector space. In this case the equation for the energy /momentum relation (29) becomes a 
matrix equation whose eigenvalues represent that energies of the neutrino energy eigenstates 
which in general are not the same as the flavor, mass and k-essence eigenstates. A formula 
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for the momentum as a function of energy can be derived from the information found in 
section §2.3. If we expand this formula around g = and m 2 = to leading order, and then 
replace g and M 2 = Q 2 m 2 with matrices that act in the neutrino flavor space we end up 
with 



P 



EI 



EG M 2 Vt~ 2 



IE 



(54) 



where G = f 



+ n ■ V0) 2 . Here G and M 2 are operators in the flavor space. In general 
it will not be possible to diagonalize G and M 2 at the same time, since there is no reason 
to assume that the k-essence and mass eigenstates are the same. We can relate the flavor, 
mass, and k-essence eigenstates by two SU(2) matrices. Define two SU{2) matrices Uf a and 
Vf m such that 



m 



Vf m \v m ), \u f ) 



Ufa Wa) 



\ u a) 



f,m 



fm 



(55) 



In this case where there are three sets of neutrino eigenstates, one can not simply disregard 
the phases in the k-essence and mass mixing matrices Uf a and Vf m as can be done in the 
standard treatment of purely mass-induced neutrino oscillations. While we can eliminate 
most of the phases through a redefinition of the neutrino eigenstates, there is not enough 
freedom to get rid of them all. After any redefinition of the neutrino eigenstates, there will 
still be a single overall phase left. Let's call this residual phase a and define the k-essence 
and mass mixing matrices as 



U 



e ia cos 9 K 
—e ia sin 9 K e ta cos 9 K 



e %a sin 9k 



V 



cos 9 m 
— sin 9 m 



sin 9 m 

cos 9 m 



If G\ and Gi are the eigenvalues of G, the momentum operator can be written as 



P 



E 



1 
1 



E. 



-U 



G x 
G 2 



17+ 



2E 



-V 



~,,2 



m 











T-,2 







(56) 



(57) 



where m 2 = f2 m~. If we substitute this for the momentum operator in the integrand in 
(38) we find that modulo terms proportional to the identity matrix we get 



EI -p 







Am 2 


— cos 29m 


sin 29m 


EAG 


- cos 29 K e~ 2ia sin 29 K 


_dt_ 





AE 


sin 29 M 


cos 29 m 


4 


e 2ia sin 29 K cos 29 K 



^2 



Am 

~Je 



— cos 29m — V cos 29k sin 29 M + y sin 29 K e 2ia 
sin 29 M + y sin 29 K e 2ta cos 29 M + y cos 29 K 



(5* 



where Am 2 

?2 



m 



mf. Here we have introduced a new parameter y, which is defined as 
y = ~ E A ^? where AG = G 2 — G\. This operator can be written in a much more compact 
form by defining new variables: 



AM 2 
AE 



i -, 2 



AM 2 
AE 



cos 29] 



sm29r 



Am 

AE 
Am 2 

AE 



(cos 29 m + y cos 29 k) 
(sin2^ M + ysin2^e" 2lQ ) 



(59) 
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Note that the phase cr, while present here, will not effect our final result. Inverting these 



AM 2 
AE 



^21 



| Am 
AE 



1 + 2?/ cos 20 + y 1 



sin 2 29 L 



sin 2 29 M + 2?/ sin 26> M sin 29 K cos 2a + y 2 sin 2 29 K 



tana 



1 + 2ycos26 + y 2 
sin 29k sin 2a 
sin 2#a/ + y sin 29k cos 2a 



where 



cos 26 = cos 29 u cos 29k + sin 2#m sin 29k cos 2a 
With these, equation (58) can be written as 



-cos2# L e - i<T sin2# L 
e ia sin 29 L cos 2fl L 





dx. 




AM 2 " 


p- 


dt_ 





AE 



(60) 



(61) 



(62) 



The phase operator is therefore 



- cos 2# L e" ifT sin 2# L 
e iCT sin 29 L cos 2# L 



cos<^(£) 



1 
1 



(63) 



where tp{£) = f* d ^fdx and £ = x d — x e . Using (63) as our neutrino evolution operator, 
the flavor eigenstates at a later time after the neutrino have traveled a distance £ are given 
by: 



\v e (t)) = isin<f(£) (co$29 L \v e ) - e m sin20 L \v^) + cos <£>(£) |i/ e ) 



M*)> 



sin </?(£) (e iCT sin26 l z / |z/ e ) + cos26>l + cosy?(£) (z/^) 



(64) 
(65) 



Therefore, if a neutrino of flavor / is created and travels a distance £ to a detector, the 
probability that the neutrino is observed as a /' flavored neutrino is 



5 



fj' 



'-iff J' sin 2 29 L sin 2 cp{£) 



(66) 



It is interesting to consider the case when the k-essence field does not vary rapidly in space 
and time compared to terrestrial scales. In that case we can treat the integrand as a constant, 
and therefore the phase becomes 



<p(£) 



^21 



| Am 
AE 



-£y / l + 2ycos2<d + y 2 



(67) 



Since the phase is proportional to the distance between the neutrino emitter and detector 
then we can write the oscillation probability as 



P(u f -> Vf >) = S fJI + (-1) 5 /./' sin 2 29 L sin : 



(61 
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where A is the oscillation length and is defined as 

tt£ A-kE 1 , , 

<p(t) I Am 2 1 v /l + 2ycos26 + ?/ 2 

It is interesting to compare results from the KINO mechanism with those from neutrino 
oscillations induced by mass. In the case where neutrino oscillations are due to mass entirely, 
the inverse oscillation length goes like 

*rA- = (70) 

Likewise, if the mass of neutrinos vanishes then the inverse oscillation length goes like 

(92 ~ 9i)u 



AnX- 1 = E\AG\ = E 



X 



n • V(f>y (71) 



Comparing (70) and (71) it is apparent that neutrino oscillations induced by either a flavor- 
non-diagonal mass term or a flavor-non-diagonal k-essence coupling will lead to noticeably 
different energy dependences for the oscillation length. If neutrino oscillations are entirely 
induced by k-essence then the oscillation length goes like A _1 ~ E. The result (71) should be 
compared to neutrino oscillations induced by either the VLI mechanism or VEP mechanism. 
Both the VLI and VEP mechanisms have the same A" 1 ~ E behavior that the KINO 
mechanism has. This should come as no surprise since the flavor-dependent emergent metric 
G$ can be viewed as a regular space-time metric but with a flavor dependent, and therefore 
equivalence principle violating, gravitational constant. Equations (68) and (71) immediately 
tell us that in order for there to be flavor oscillations not only must g\ ^ g 2 , but also u ^ 0. 
Since u vanishes for c s = 1 this means that there must be lorentz violation if k-essence is 
to have any effect on neutrino oscillations. K-essence, therefore, acts as a lorentz violating 
"aether" background, and can be the motivation behind models of neutrino oscillation that 
invoke Lorentz violation. 

Analysis of the available data from current and past neutrino observatories have tended 
to favor mass induced neutrino oscillations, which can produce the desired A oc E type 
behavior. Therefore, the KINO mechanism alone does not suffice to explain the observations 
seen in the numerous neutrino experiments that have been carried out. However, it is still 
possible that k-essence could be a subleading contribution to neutrino oscillations with mass 
being the dominant cause. Studies have looked into the possibility of alterations to the 
leading order A oc E dependence and have been able to place very tight constraints on the 
coefficients of sub-leading contribution to the energy dependence of A. In [11] they considered 
the possibility of different mechanisms inducing oscillations in the <-> v T channel, among 
which were the VLI and VEP mechanisms. Since the VLI and VEP scenarios both lead to 
the same energy dependence for the oscillation length that KINO does, their constraints on 
the VLI and VEP coefficients can be easily translated to a bound on AG 

|AG| < 1.2 x 10~ 23 ( 72) 

this is the most conservative bound that can be placed on |AG| and is independent of the 
mixing angle 6 k- This bound would seem to cast doubt on KINO as an even subleading 



18 



effect in the v yi <-» v T channel. We would like to emphasize that this is a strong indication of 
the nature of the coupling of neutrinos to the k-essence scalar. Our analysis suggests that if 
the k-essence scalar field exists, in order to be phenomenologically viable, its couplings must 
be flavor blind. 

Although KINO would seem to be immediately discounted from consideration, it may be 
possible to realistically consider this model if we are ready to include further symmetry 
violating terms in our action. It has been shown in [12] that if lorentz and certain types of 
CPT violating terms are included in the neutrino action then it becomes possible to create 
pseudo-mass terms at high energies just when lorentz violation effects should be taking over. 
These types of models open up the possibility of a unified explanation for all the existing 
neutrino data, including the controversial LSND results [43]. Since these models have been at 
least qualitatively compatible with experiment, this leaves open the possibility that k-essence 
could still play a role in neutrino oscillations. 

As we have just mentioned this approach requires a specific kind of CPT violating term in 
the neutrino action. Although it is always possible to put terms into the action arbitrarily, 
in principle it might be possible for k-essence to be a source for these as well. K-essence 
could lead to the needed CPT breaking terms in the neutrino action by one of two ways, the 
first is by a possible axial vector term in the spinor covariant derivative 6 . In this paper we 
have assumed that space-time is flat, and as a result the covariant derivative is proportional 
to 7^. However, if more general curved space-times are considered then in general the spinor 
connection will have a non-zero axial vector part. Another way in which k-essence can 
source the necessary CPT violating term is by considering the possibility of torsion in the 
emergent space-time. Such a treatment of k-essence would require taking into account of the 
spin-orbit coupling between the neutrino and the k-essence field and treating the metric and 
connection on the emergent space-time as independent variables each requiring their own 
field equations. Torsion in k-essence backgrounds is an interesting possibility and deserves 
further study in its own right. 

4 Conclusion 

In this paper we have investigated the effects on neutrino velocity and oscillations, when neu- 
trinos couple to a k-essence field background though the emergent metric G^ u . Specifically 
we have studied the results of replacing the vierbein of the gravitational background with 
the vierbein E% of the emergent k-essence geometry in the neutrino action. The first impli- 
cation of this coupling is a change in the neutrino dispersion relation, which means that the 
neutrino velocity would be dependent on the k-essence background in which it propagates. 
Without k-essence massless neutrinos always propagate at the speed of light with respect 
to local observers. However, with k-essence, massless neutrinos will move at a new speed 
c v with respect to an observer who perceives a uniform k-essence background. Therefore, if 
future observations show that neutrinos do travel at less than the speed of light, it cannot be 
conclusively determined by this alone if this is due to neutrinos being massive or if neutrinos 

6 See appendix A 
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are massless but coupled to k-essence. However, at present, no observations have found a 
measurable difference between the speed of light and neutrinos. If neutrinos are assumed to 
be massless, than past observations constrain the deviation of c v from the speed of light to 
no more than 1 part in 10 8 . 

The other effect that k-essence has on neutrinos is in the phenomenon of neutrino oscilla- 
tions. We have shown that if the k-essence-neutrino coupling g is promoted to the status of 
an operator g that acts in the neutrino flavor space, neutrino oscillations are produced even 
in the absence of a neutrino mass term. We found that if neutrino oscillations are caused 
entirely by k-essence, then the oscillation length depends on the energy like A ~ E~ x . This is 
to be contrasted with the result from mass induced neutrino oscillations where A ~ E. Thus, 
while neutrino oscillations can be induced by k-essence, it will lead to a drastically different 
observation of the energy dependence of the oscillation length. Current data seems to favor 
aA~£ behavior which implies an important constraint on the couplings of neutrinos to the 
k-essence scalar. In order to be phenomenologically viable these couplings must be flavor 
diagonal. Our analysis is a very good example of how interactions between dark energy and 
visible matter can be used to constrain their couplings. 

Neutrino oscillations induced by k-essence have many of the same properties that some 
[19, 20, 18] earlier proposed mechanisms had, and in fact k-essence can be seen as a realization 
of these past phenomenologically motivated models. In VLI models the lorentz violation is 
incorporated by introducing some preferred 4-vector into the neutrino action. This preferred 
4- vector has the interpretation in k-essence as the 4-gradient of the k-essence field. In fact, 
as can be seen from the formulae in section §3, the physical effects vanish at c s = 1 where 
our model has exact Lorentz invariance. The VEP mechanism attempts to explain the 
source of neutrino oscillations as a consequence of a flavor- non- diagonal coupling to gravity. 
Although perhaps correct from a theoretical point of view, this mechanism calls for a very 
drastic change in our understanding of fundamental physics; namely it requires us to give 
up the long held notion of the equivalence principle. The type of coupling studied in this 
paper mimics exactly the VEP mechanism, but since the flavor non-diagonal coupling is 
in the k-essence sector instead of the the gravitational sector, the equivalence principle is 
maintained. In short, k-essence can be seen as the source of lorentz violation in the VLI 
model, and alternatively, a reinterpretation of the VEP model. KINO can therefore be seen 
as a theoretically and phenomenologically motivated manifestation of earlier attempts at 
alternatives to mass-induced neutrinos oscillations. 

In this paper, we have coupled the neutrinos to the k-essence scalar using the emergent 
metric which is covariantly constant. In a future work we will show how this can be gen- 
eralized to more general couplings by introducing the analogue of torsion for the emergent 
space-time. This generalized coupling will allow us to analyze the role, if any, of CPT vio- 
lation in the coupling of k-essence to fermions. If CPT violation can be naturally produced, 
then it would be possible to discuss more realistic models of neutrino oscillations as discussed 
in [42]. 
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A Dirac Fermions in the Emergent K-essence Metric 

The simplest action of a Dirac fermion ip coupled to a metric is conventionally written 

as 



S = J d A xEi) [ij a E^V^ -M\if> 



(73) 



where E% is the vierbein corresponding to the metric G^, and is the spinor covariant 
derivative which is given by 



T^n — dfj, — -£l a bnO~ ab . 





4' 



Here a ab = | [7 a ,7 6 ] and fl a b^ is the spinor connection form which by definition is 



ab/i 



Eaud^El + E av E°T 



(74) 



(75) 



where Y v a is the standard Christoffel symbol. After some work, the spinor covariant deriva- 
tive can be shown to satisfy 



where 



1 



d, + \El (d,E b v - d v El) + Z - l5 A, 



A^ = -e abcd E afl (d a E bu — d u E ba ) E v c E c d 



(76) 



(77) 



The formula (76) is valid not just for the emergent geometry of k-essence, but for all geome- 
tries with a vierbein. In this paper we will assume that the space-time geometry is flat. In 
this case the axial vector part vanishes since = always if space-time is flat, and the only 
non-zero part of the spin connection is the vector portion. Therefore, the spinor covariant 
derivative for a general k-essence field and a flat space-time background is 



Tp = 



gd u u 



2(1 + gu) 



+ 9U a x (ft dU(pd ^ 



3 d„tt f gu 
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It the cosmologically relevant case where is time dependent but Vj0 = 0, the covariant 
derivative in this case leads to 



In k-essence models that attempt to explain the cosmological constant problem, the higher 
order derivatives of the field become irrelevant at late times since the field reaches a steady 
state by that point. Therefore, it is reasonable to ignore the spinor connection in our analysis. 
However, if one is interested in the effect k-essence has on neutrinos in the early universe, in 
particular around the time of matter-radiation equality, then the spinor connection can be 
an important contribution to the neutrino action. 
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